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In the cosmological paradigm, cold Dark Matter (DM) dominates the mass content of the Universe 
and is present at every scale. Candidates for DM include many extensions of the standard model, 
with a weakly interacting massive particle (WIMP) in the mass range from 50 GeV to greater than 
10 TeV. The self-annihilation of WIMPs in astrophysical regions of high DM density can produce 
secondary particles, including very high energy (VHE) 7-rays, with energies up to the dark matter 
particle mass. The VERITAS array of Cherenkov telescopes, designed for the detection of VHE 
7-rays in the 100 GeV- 10 TeV energy range, is an appropriate instrument for the detection of DM 
and is complementary to Fermi-LAT. Dwarf spheroidal galaxies (dSphs) of the Local Group are 
potentially the best targets to search for the annihilation signature of DM due to their proximity 
and large DM content. We report on the latest VERITAS observations of dSphs and discuss the 
results in the framework of WIMP models. 



I. INTRODUCTION 

The compelling evidence for the presence of DM in the 
different structures of the Universe [l[ has motivated 
numerous efforts to search for DM by means of as- 
trophysical observations. If DM is made of WIMPs 
annihilating into standard model particles, indirect 
searches for DM annihilations with VHE 7-rays pro- 
vide a very promising way to constrain the nature of 
the DM particle: 7-rays are free of any propagation 
effects on short distances (< IMpc) and DM parti- 
cle annihilation is predicted to give a unique 7-ray 
spectrum. Such searches are often conducted using 
pointed observations toward nearby DM overdensi- 
ties, because the annihilation rate is proportional to 
the squared DM density. Popular targets include the 
Galactic Center |2|-t5jl , satellite galaxies of the Milky- 
Way (MW) [i^Sl globular clusters [ll[l3 and clus- 
ters of galaxies 

The dSphs of the Local Group best meet the crite- 
ria for a clear and unambiguous detection of DM. 
They are gravitationally bound objects and contain 
up to O(IO^) more mass in DM than in visible mat- 
ter. As opposed to the Galactic Center, they are en- 
vironments with a favorably low astrophysical 7-ray 
background. Neither astrophysical 7-ray sources (su- 
pernova remnants, pulsar wind nebulae, etc) nor gas 
acting as target material for cosmic rays have been 
observed in these systems [l6l | . Furthermore, their rel- 
ative proximity and high galactic latitude make them 
the best astrophysical targets for high signal-to-noise 
detection. 

This paper gives a status of the dSph observational 
program carried out by the ground-based VHE 7-ray 
observatory VERITAS. Section |lT] presents the instru- 
ment and summarizes the VERITAS dSph observa- 
tions and data analysis. Section IlIII interprets the re- 



sults in terms of constraints on DM models, focusing 
on models which include a Sommerfeld enhancement. 
Finally, section ITVl is devoted to the conclusion. 



II. VERITAS OBSERVATION AND 
ANALYSIS 

VERITAS is an array of four 12-meter imaging atmo- 
spheric Cherenkov telescopes (lACTs) located at the 
base camp of the F. L. Whipple Observatory in south- 
ern Arizona. Each VERITAS telescope consists of a 
large optical reflector which focuses the Cherenkov 
light emitted by particle air showers onto a camera 
of 499 photomultiplier tubes. The large effective area 
(~ 10^ nr^), in conjunction with the stereoscopic imag- 
ing of air showers, enables VERITAS to be sensitive 
over a wide range of energies (from 100 GeV to 30 
TeV) with an energy and angular resolution of 15- 
20% and 0.1° per reconstructed 7-ray, respectively. 
VERITAS is able to detect a point source with 1% 
of the Crab Nebula flux at a statistical significance 
of 5 standard deviations above background (5 cr) in 
approximately 25 hours of observations. For further 
details about VERITAS, see, e.g., [l3|. 
Since 2007, the VERITAS observatory has started an 
active observational program on a dSph sample rep- 
resenting the most promising targets visible in the 
northern hemisphere. Table |T] presents the current 
observations carried out by VERITAS on five well- 
studied and closeby dSphs, namely the Draco, Ursa 
Minor, Bootes 1, Willman 1 and Segue 1 galaxies. 
The so-called wobble pointing mode, where the point- 
ing position is shifted by ± 0.5 deg from the target 
position, has been used during these observations. 
The wobble pointing strategy allows for simultane- 
ous background estimation and source observation. 
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dSph target 


Draco 


Ursa Minor 


Bootes 1 


Willman 1 


Segue 1 


Exposure [hours] 


18.4 


18.9 


14.3 


13.7 


51.5 


Zenith angle [deg] 


26-51 


35-46 


17-29 


19-28 


16-39 


Excess [counts] 


-28.4 


-30.4 


28.5 


-1.45 


30.4 


Significance [a] 


-1.5 


-1.8 


1.3 


-0.1 


0.9 


95% CL upper Umit [counts] 


18.8 


15.6 


72.0 


36.7 


135.9 


^95%CL(g > 300 GeV) [cm-'s-i] 


0.5 X 10"^^ 


0.4 X 10~^^ 


2.2 X 10"^^ 


1.2 X 10"^^ 


0.8 X 10"^^ 



TABLE I: The VERITAS dSph observations and corresponding analysis results. The significance is calculated according 
to the method of Li & Ma [13] . The 95% CL ULs on the number of 7-rays in the ON-source region is derived using the 
Rolke prescription The ULs on the integrated flux have been computed assuming a power law of index -2.6. 



thus reducing the systematic uncertainties in the back- 
ground determination. Data reduction follows the 
methods described in 18|. After calibration of the 



sight (s) and over the solid angle Af2: 



data, the 7-rays are selected by calculating the Hillas 
parameters of the recorded camera images. The se- 
lection cuts were optimized for the detection of a 5% 
Crab Nebula-like source. After the 7-ray selection, 
the residual hadronic background in a circular region 
of 0.12 deg radius centered on the target position 
(called the ON-source region) is estimated using the 
reflected background model. The analysis of the data 
did not reveal any significant excess over the estimated 
background in any of these observations. Table U dis- 
plays the measured 7-ray excesses for each of the dSph 
dataset, along with the corresponding significances, 
the resulting upper limits (ULs) on the number of 7- 
rays in the ON source region and on the integrated 
flux above 300 GeV. The 95% confidence level (CL) 
integral flux ULs lie in the range 0.4-1.2% of the Crab 
Nebula flux. 



III. CONSTRAINTS ON DARK MATTER 
MODELS 

A. Classical constraints 

The absence of signal in any of these observations can 
be used to derive constraints on various dark matter 
models. The 7-ray flux from the annihilations of DM 
particles, of mass moM, in a spherical DM halo is 
given by a particle physics term times an astrophysics 
term: 



dE 



(Arj,E) 



(crv) dN^ 



dE 



xj(Ar2). (1) 



The particle physics term contains all the informa- 
tion about the DM particle: its mass m^Mj its to- 
tal velocity-weighted annihilation cross-section (crv} 
and the differential 7-ray spectrum from all flnal 
states weighted by their corresponding branching ra- 
tios, dN^/dE. The astrophysical factor J(Ar2) is the 
square of the DM density integrated along the line of 



J(A^I) = / 

J An 



p2(r[s])ds. 



(2) 



los 



The solid angle is given here by the size of the signal 
search region defined previously in the analysis, i.e. 
6 < 0.12 deg. The estimate of the astrophysical fac- 
tor requires a modeling of the dSph DM distribution. 
Each dSph DM distribution has been modeled with a 
Navarro, Frenk & White (NEW) profile except 
Segue 1 for which the DM distribution has been mod- 
eled with an Einasto profile [22^ . The parameters of a 
dSph profile are constrained using its star kinematics 
and have each been taken from the literature 0, . 
The resulting astrophysical factors of Draco, Ursa Mi- 
nor, Bootes 1, Willman 1 and Segue 1 are presented 
in Table ini along with their distance and profile mod- 
eling. 

Once the astrophysical factors have been estimated, 
the ULs in the number of 7-ray candidates for each 
dSph can be translated into ULs on the total annihi- 
lation cross-section: 



(av) 



95% CL 



Stt 



DM 



Tob./o™Aff(E)^|^dE' 



(3) 

where Tobs is the total observation time and ^eff (E) is 
the effective area as a function of energy, zenith angle 
and the offset of the source from the pointing posi- 
tion. 

Figure n shows the VERITAS 95% CL ULs on the 
velocity-weitghted annihilation cross-section (crv) as 
a function of the DM particle mass for each of the 
five observed dSph listed in Table HI The Draco, Ursa 
Minor, Bootes 1 and Willman 1 limits have been ob- 
tained using a composite DM spectrum with a 90 % 
bb branching ratio and a 10% t+t~ branching ratio 
0] ■ The Segue 1 limits have been obtained assuming a 
pure W"*'W~ DM annihilation spectrum, which is sim- 
ilar to the composite bb - t'^t~ spectrum. The limits 
are in the range 10^^^ — 10"^^ cm'^ s~^, depending on 
the considered dSph. The best limits are obtained 
for the Segue 1 dSph, which has the deepest exposure 
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dSph target 


Draco 


Ursa Minor 


Bootes 1 


Willman 1 


Segue 1 


Profile modeling 


NFW 


NFW 


NFW 


NFW 


Einasto 


Distance [kpc] 


80 


66 


62 


38 


23 


J(6l < 0.12°) [GeV^ cm~^ sr] 


1.5 X 10^* 


2.7 X 10" 


1.1 X 10" 


8.4 X 10" 


7.7 X 10" 



TABLE II: Profile modeling, distance and astrophysical factors of tlie five dSphs observed witli VERITAS. 




boost factors up to 10^. In such scenarios, the anni- 
hilation cross-section is given by: 



10" 

[GeV] 



FIG. 1: 95% CL ULs on (crv) as a function of the DM 
particle mass for the five dSphs observed with VERITAS . 
The limits are computed using a composite DM annihila- 
tion spectrum (10% t^t" + 90% bb), except for the case 
of Segue 1 in which a pure W""" W^ channel was considered. 
The grey band area represents a range of generic values for 
the annihilation cross-section for thermally produced DM. 



and one of the highest astrophysical factor among the 
five dSphs (see table . The 95% CL upper limit on 
(ctv) is at the level of 8 x 10^^'' cm'^ s~^ around 1 TeV 
and is among the best reported so far with dSph VHE 
observations. 



B. Constraints on models with a Sommerfeld 
enhancement 

In this section, we consider DM models where the an- 
nihilation cross-section can be boosted with the Som- 
merfeld enhancement effect. We use the Segue 1 data 
to constrain such models. The Sommerfeld enhance- 
ment is a particle physics effect that was suggested to 
explain the cosmic ray lepton anomalies measured by 
ATIC and PAMELA [2i,[2i|. It arises when the two 
DM particles interact through an attractive potential, 
mediated by the exchange of a massive boson 
The annihilation cross-section enhancement is partic- 
ularly effective when the DM particle relative velocity 
is very small. Depending on the mass m^ and the 
coupling a of the exchanged boson, the Sommerfeld 
enhancement can exhibit a series of resonances for spe- 
cific values of the DM particle mass, giving very large 



(crv) = S X (crv)o, 



(4) 



where S is the Sommerfeld boost factor averaged over 
the DM particle velocity distribution and depends on 
the DM particle mass, (o'v)o is the thermal WIMP 
annihilation cross-section before freeze-out. The Som- 
merfeld enhancement is of particular interest for cold 
DM halos like dSphs, where the mean star velocity 
dispersion can reach a few kms""'^. 
We focus here on two models comprising a Sommer- 
feld enhancement to the annihilation cross-section. 
The first model (hereafter model I, [131) assumes that 
the dark matter particle is a wino-like neutralino 
(arising in SUSY extensions of the standard model), 
which annihilates through the exchange of a TP boson 
(mzo ^ 90 GeV, a ^ 1/30), thus leading to a Som- 
merfeld enhancement. The second model (hereafter 
model II) introduces a new force in the dark sector 
psj . carried by a light scalar field (j) predominantly 
decaying into leptons and with a mass 0(1 GeV). In 
such models, dark matter annihilates to a pair of 
(f> scalar particles, with an annihilation cross-section 
boosted by the Sommerfeld enhancement. The cou- 
pling a of the light scalar particle (j) to the dark mat- 
ter particle is determined assuming that XX ~^ 4"^ 
is the only channel that regulates the dark matter 
density before freeze-out [2^. Figure shows the 
VERITAS constraints on each of these models, de- 
rived with the observations of Segue 1. The left 
panel of Figure [5] shows the constraints on model 
I. The Sommerfeld enhancement exhibits two reso- 
nances in the considered dark matter particle mass 



for 



4.5 TeV and m^ ~ 17 TeV, respec- 



range, 

lively. VERITAS excludes these resonances, which 
boost the annihilation cross-section far beyond the 
canonical (crv) 3 x 10~^^ cm"^ s~^. The right panel 
of Figure [2] shows the VERITAS constraints on model 
II, for a scalar particle with mass m^ = 250 MeV. Two 
channels in which the scalar particle decays either to 
e"''e~ or fi^fJ," have been considered. VERITAS ob- 
servations start to disfavor such models, especially 



for the e"'"e" 

nances are beyond (crv) 3 x 10~^^ cm'^ s~^. This re- 
sult holds for (p particle masses up to a few GeV. 



channel where some of the reso- 
"^^ cm^ s" 
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FIG. 2: 95% CL exclusion curves from the VERITAS observations of Segue 1 on ((tv)/S as a function of the dark matter 
particle mass, in the framework of two models with a Sommerfeld enhancement. The expected Sommerfeld enhancement 
S applied to the particular case of Segue 1 has been computed assuming a Maxwellian dark matter relative velocity 
distribution [3] . The grey band area represents a range of generic values for the annihilation cross-section in the case of 
thermally produced dark matter. Left: model I with wino-like neutralino dark matter annihilating to a pair of W'''W~ 
bosons. Right; model II with a 250 MeV scalar particle decaying into either e+e~ or /x"'"/^". 



IV. SUMMARY AND CONCLUSION 

Since its commissioning phase, the VERITAS obser- 
vatory has been actively observing the best dSph 
targets visible from the northern hemisphere, with 
data currently taken toward the Draco, Ursa Mi- 
nor, Bootes 1, Willman 1 and Segue 1 galaxies. 
The non-detection of a VHE 7-ray signal in any of 
these observations has lead to limits on the velocity- 
weighted annihilation cross-section (crv) of the order 
of 10~^^ — 10~^^ cm^ s~^. The Segue 1 data have pro- 
vided the best limits on (crv) reported so far with 
any dSph observations, though these limits are at 
least two orders of magnitude away from the canoni- 
cal (ctv) ~ 3 X 10~^^cm3s~^. The Segue 1 data have 
also been used to constrain models with a Sommerfeld 



enhancement, disfavoring them as being a good inter- 
pretation of the cosmic ray lepton excesses measured 
by ATIC and PAMELA. 
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